Evolution of the 21 cm signal throughout cosmic history 
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The potential use of the redshifted 21 cm line from neutral hydrogen for probing the epoch 
of reionization is motivating the construction of several low-frequency interferometers. There is 
also much interest in the possibility of constraining the initial conditions from inflation and the 
nature of the dark matter and dark energy by probing the power-spectrum of density perturbations 
in three dimensions and on smaller scales than probed by the microwave background anisotropies. 
Theoretical understanding of the 21 cm signal has been fragmented into different regimes of physical 
interest. In this paper, we make the first attempt to describe the full redshift evolution of the 21 
cm signal between < z < 300. 

We include contributions to the 21 cm signal from fluctuations in the gas density, temperature 
and neutral fraction, as well as the Lya flux, and allow for a post-reionization signal from damped 
Lya systems. Our comprehensive analysis provides a useful foundation for optimizing the design of 
future arrays whose goal is to separate the particle physics from the astrophysics, either by probing 
the peculiar velocity distortion of the 21 cm power spectrum, or by extending the 21 cm horizon to 
z > 25 before the first galaxies had formed, or to z < 6 when the residual pockets of hydrogen trace 
large scale structure. 

PACS numbers: 98.80 Hw, 95.35.+d, 98.62.-g 
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I. INTRODUCTION 

Modern cosmology has advanced considerably in re- 
cent years due to precise measurements of the cosmic 
microwave background (CMB) and large-scale structure 
[l[ . These data sets probe the distribution of matter at 
redshifts z ~ 10 3 and z < 0.3, respectively, and only 
cover ~ 0.1% of the entire comoving volume of the ob- 
servable Universe Q. In addition, significant informa- 
tion about individual galaxies and quasars out to z = 7 
has been obtained through the efforts of many telescopes 
at multiple wavelengths Q. Nonetheless, our knowledge 
and understanding of the evolution of most of the bary- 
onic matter, which resided in the intergalactic medium 
(IGM), is limited. Constraints from WMAP3 Q on the 
optical depth to the CMB (r = 0.093±0.029) only weakly 
constrain the redshift of reionization. The presence of a 
Gunn-Peterson trough [j| in quasar spectra at z > 6.5 
provides a loose limit on the end of reionization Q . Ab- 
sorption by the Lya forest in quasar spectra constrain 
the IGM temperature to be T K = 2 x 10 4 K at z ~ 4 Q. 
Beyond these fragmentary data points lies a complicated 
story of ionization and heating of the IGM going back to 
the formation of the first stars [1, Q . 

Low frequency observations of the redshifted 21 cm line 
of neutral hydrogen present one of the most promising av- 
enues for exploring the entirety of the poorly constrained 
period between recombination and reionization [lol [ll[ . 
Additionally, mapping hydrogen throughout this period 
and at lower redshifts raises the prospect of measuring 
primordial density fluctuations over a much larger vol- 
ume than that probed by the CMB or galaxy redshift 



surveys [12j, [l3(. Several experiments are currently be- 
ing constructed (such as MWA 1 , LOFAR 2 , PAPER 3 , 
21CMA 4 ) and more ambitious designs are being planned 
(SKA 5 ) to detect the theoretically-predicted 21 cm sig- 
nal. 

The volume to be mapped can be loosely divided into 
three parts: (i) the "Dark Ages" between the epoch of 
cosmic recombination and the appearance of the first 
galaxies (25 < z < 10 3 ); (ii) the epoch of reionization 
(EoR) during which hydrogen was ionized 6 by the UV 
radiation from galaxies (6 < z < 25); and (Hi) the post- 
reionization epoch (z < 6). The transition between the 
first two epochs, characterized by the appearance of the 
first luminous sources, is a period of "Cosmic Twilight'' , 
during which the radiation fields of the sources begins to 
affect the gas, heating and exciting the largely neutral 
inter-galactic medium (IGM). Eventually, the small ion- 
ized regions surrounding groups of sources expand and 
percolate beginning the process of reionization. After 
reionization, the residual neutral hydrogen is found in 
dense clumps (such as damped Lya systems [15J), self- 
shielded from ionizing radiation. 

Precise measurements of the power-spectrum of 21 cm 
brightness fluctuations could constrain the initial condi- 
tions from inflation fl2| (including deviations from scale- 
invariance or Gaussianity) , as well as the nature of the 
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However, exotic physics such as the decay of dark matter [13 ] 
or evaporation of primordial black holes might alter this picture 
dramatically. 
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dark matter (including a measurement of the neutrino 
mass in the range expected from atmospheric neutrino 
experiments 0|) and dark energy [lq ]. The 21 cm fluctu- 
ations are expected to simply trace the primordial power- 
spectrum of matter density perturbations either before 
the first galaxies had formed (z > 25) [HI, H3] or after 
reionization (z < 6) - when only dense pockets of self- 
shielded hydrogen survive [H, [l9[ . During the epoch of 
reionization, the fluctuations are mainly shaped by the 
topology of ionized regions [13, HH [HJ , and thus depend 
on astrophysical details. Nevertheless, even during this 
epoch, the line-of-sight anisotropy of the 21 cm power 
spectrum due to peculiar velocities, can in principle be 
used to separate the implications for fundamental physics 
from the unknown details of the astrophysics [13, H3| • 

Recent work has explored the 21 cm signal that may 
arise during reionization [24], HE} from variations in the 
Lya flux H3, X-rayheating [H[ , gas density 
fl2l H3] or peculiar velocities [20l I23L I30I ] . In addition, the 
post-reionization signal has been explored by Wyithe and 
Loeb [18j. In general, previous studies have focused on 
individual regions of parameter space where only one of 
two of these fluctuations needs to be considered and the 
rest are assumed to be negligible. Here we attempt for 
the first time to present a global picture of the redshift 
evolution of the 21 cm signal, including all relevant fluc- 
tuations. The full history can be used to select particu- 
lar redshift (frequency) windows and observing strategies 
within which fundamental physics would be optimally 
constrained. 

In analyzing the full history of 21 cm fluctuations, we 
seek to address two main issues. First, we would like 
to identify the redshift at which the presence of astro- 
physical sources begin to modify the pristine cosmologi- 
cal signal. This defines the minimum redshift that must 
be accessible for a future observatory aiming to detect 
directly the pristine signal (in analogy with CMB exper- 
iments). Although the answer is model-dependent, we 
will show that the rapidity with which structures form 
once non-linear collapse gets underway means that the 
model dependence is weak. The weak dependence we find 
is generic to a Gaussian random field of initial density 
perturbations, as expected from inflation, for which the 
fraction of mass incorporated in collapsed objects grows 
exponentially (on the high-density tail of the Gaussian) 
at early times. Second, we wish to address the possibility 
of using the angular separation proposed by Barkana and 
Loeb [23| to recover cosmological information from 21 cm 
observations after astrophysical effects become large. We 
will show that there exists a window in redshift where 
this technique might be used to extract the dark-matter 
power spectrum despite the complication of other sources 
of fluctuations. Finally, we would like to identify the pri- 
mary astrophysical parameters that control the 21 cm 
signal at subsequent epochs. 

The paper is organized as follows. In ijTT]we discuss the 
calculation of the 21 cm signal, attempting to bring to- 
gether a variety of different sources of fluctuations. The 



results of these calculations are illustrated in ^IIII In 
QTVl we discuss the experimental requirements for de- 
tecting the 21 cm signal and explore the performance of 
three fiducial experiments. Finally, we conclude in ifVl 
and discuss future prospects. 

Throughout this paper, we assume a cosmology with 
n m = 0.26, Q A = 0.74, fl b = 0.044, H = 
lOO/ikms^Mpc" 1 with h = 0.74, n s = 0.95, and 
erg = 0.8, consistent with the latest measurements Q. 



II. THEORETICAL FRAMEWORK 

A. 21 cm brightness temperature 

We begin by briefly summarizing the physics of the 21 
cm signal and refer the interested reader to Furlanetto, 
Oh, and Briggs for further details. The 21 cm line of 
the hydrogen atom results from hyperfine splitting of the 
IS" ground state due to the interaction of the magnetic 
moments of the proton and the electron. The HI spin 
temperature Ts is defined through the ratio between the 
number densities of hydrogen atoms in the IS triplet and 
IS singlet levels, ni/no = (gi/go)exp(—T is /Ts), where 
{91/90) = 3 is the ratio of the spin degeneracy factors 
of the two levels, and T* = hcjk\i\ cm = 0.068 K. The 
optical depth of this transition is small at all relevant 
redshifts, yielding brightness temperature fluctuations 
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Here Xm is the neutral fraction of hydrogen, S b is the 
fractional overdensity in baryons, and the factor of 4/3 
arises from including the effect of peculiar velocities and 
averaging the signal over angles. The spin temperature 
is given by 



rp-l _ 1 + X a T a 1 + X C T K 

1 "f" 3>Q ~\~ Xq 



(2) 



where T a is the color temperature of the Lya radiation 
field at the Lya frequency and is closely coupled to Tk by 
recoil during repeated scattering, and the cou- 

pling coefficients to atomic collisions and Lya photons, 
respectively. The spin temperature becomes strongly 
coupled to the gas temperature when xtot = x c + x a > 1 . 
The collisional coupling coefficient is given by 



47; 
3AiqT 7 
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(3) 



where A\q — 2.85 x 10 -15 s _1 is the spontaneous emission 
coefficient, is tabulated as a function of T k [3l|, [32j 
and Ki? is taken from Ref. (33|. For a more detailed 
analysis of the collisional coupling, see Ref. [34T |. 
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The Wouthysen-Field effect [35|, |36j coupling is given 



_ l6Tr 2 T*e 2 f a 
27 AiqI jin e c 



(4) 



where f a = 0.4162 is the oscillator strength of the Lya 
transition. S a is a correction factor of order unity, 
which describes the detailed structure of the photon 
distribution in the neighborhood of the Lya resonance 
[37l [H, HH, H(| • We make use of the approximation for 
S a outlined in Ref. [4(|. For the models considered in 
this paper, Lya coupling dominates over collisional cou- 
pling at redshifts z < 25, although we consider this more 
carefully later. 

Fluctuations in the 21 cm signal may be expanded to 
linear order [l(| 



£>T b = Pb^b + PxSx + PaS a + Pt$T ~ S v , 



(5) 



where each Si describes the fractional variation in the 
quantity i and we include fluctuations in the baryon den- 
sity (b), neutral fraction (x), Lya coupling coefficient (a), 
gas temperature (T), and line-of-sight peculiar velocity 
gradient (v). The expansion coefficients are given by 
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It is important to realise that since fluctuations in xh can 
be of order unity, an expansion to second order becomes 
necessary near to the end of reionization. This leads to 
quartic terms contributing to the power spectrum. 

Noting that in Fourier space Sg v = —p?S [3Cj . where \x 
is the angle between the line of sight and the wavevector 
k of the Fourier mode, we may use equation (J5J to form 
the power spectrum 



Pbb + Pxx + P aa + Ptt + 2P bx 
+2P ba + 2P bT + 2P xa + 2P X T ~\ 
+PxSx8 + other quartic terms 
2p? (P u + P xS + P aS + P TS ) 
^Pss 

2Px65vX ~t~ PxS^S^x 

+other quartic terms with 5 V . 



2P aT 



(7) 



Here we note that all quartic terms must be quadratic 
in xh and separate them depending upon whether they 
contain powers of 8 V or not. Those that contain powers of 



5 V will not be isotropic and so contribute to the angular 
dependence of PT b (see [2FJ] for further discussion). 
We may rewrite Eq. ([7]) in more compact form 

P Tb (fc, /i) =P fl o(k)+ t l 2 P^(k) + p%t (k)+P mil) (k,fj). 

(8) 

where we have grouped those quartic terms with anoma- 
lous /i dependence into the term Pf(k,^)(k, fJi)- In princi- 
ple, high precision measurements of the 3D power spec- 
trum will allow the separation of Px b (fc, fi) into these four 
terms by their angular dependence on powers of /x 2 [23| . 
The contribution of the Pf(k, f i)(k, /i) term, with its more 
complicated angular dependence, threatens this decom- 
position (20l | . Since this term is only important during the 
final stages of reionization, we will ignore it in this paper 
noting only that the angular decomposition by powers of 
/m 2 may not be possible when ionization fluctuations are 
important. 

It is unclear whether the first generation of 21 cm ex- 
periments will be able to achieve the high signal-to-noise 
required for this separation [2(|. Instead, they might 
measure the angle averaged quantity 



P Tb (fc) = P„ (fc) + /> (fc)/3 + P M 4 (fc)/5, 



(9) 



(where we neglect the /x) term). In presenting 

our results, we will concentrate on Px b (fc), which is easiest 
to observe and discuss the angular separation separately 
in mil CI We will typically plot the power per logarithmic 
interval A = [fc 3 P(fc)/2 7 r 2 ] 1 / 2 . 



B. Astrophysical modeling 

Before calculating the evolution of the fluctuations, we 
must first provide an astrophysical model for the sources 
of radiation that modify the 21 cm signal. We follow the 
model of Furlanetto [4l| with some minor modifications. 

We must specify a luminosity and spectrum for sources 
of ionizing photons, Lyman series photons, and X-rays. 
The variation in spectra has relatively little effect on the 
fluctuations (but see Refs. [13] and [29[ for further discus- 
sion on this point) , so that for our purposes we will fix the 
spectrum and consider models with different luminosi- 
ties. We specify three key parameters for our astrophys- 
ical sources: the number of ionizing photons produced 
per baryon in stars that contribute to ionizing the IGM 
NionJGM = /csc^ion, the number of Lya photons pro- 
duced per baryon in stars N a , and the energy in X-rays 
produced per baryon in stars ex ■ With this parametriza- 
tion, we are separating the physics of the sources from 
the star formation history, as the quantity of physical 
relevance will be the above numbers multiplied by the 
efficiency of star formation /*. For example, the conven- 
tional definition of the ionizing efficiency £ = N- lont iGuf*i 
gives the important quantity when tracking reionization. 

For convenience, we further define f a and fx via 
N a = f a N atie { and ex = fx^x.rd where we take the 
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reference values appropriate for normal (so-called, Pop- 
ulation IP) stars N axe f — 6590 and for starburst galax- 
ies with a power-law spectrum of index as = 1.5 and 
£x,ref = 560 eV (42[. We use the spectra appropriate for 
these sources throughout. For comparison, in this nota- 
tion, the very massive (Population III) stars have [43[, 
N a = 3030 (f a = 0.46), when the contribution from 
higher Lyman series photons is included. We expect the 
value of f a to be close to unity. In contrast, fx is rela- 
tively unconstrained with values fx ^ 10 3 possible with- 
out violating the WMAP CMB optical depth constraints. 
Constraining this parameter will mark a step forward in 
our understanding of the thermal history of the IGM and 
the population of X-ray sources at high redshifts. 

Once reionization is complete we assume that photo- 
ionization heats the gas to 3 x 10 K. This is needed only 
for determining the thermal cutoff scale for the baryonic 
fluctuations for the post-reionization 21 cm signal. 



C. Calculation of fluctuations 

The focus of this paper is to combine different sources 
of 21 cm fluctuation into a complete picture of the evolu- 
tion of the 21 cm signal from the epoch when the cosmic 
gas thermally decoupled from the CMB (z ~ 200) to 
the present time. In calculating the various sources of 
fluctuations, we draw heavily upon the existing litera- 
ture. Below we briefly discuss the relevant calculations 
and then describe how we piece them together. 

The 21 cm signal has been most extensively stud- 
ied within the epoch of reionization (EoR), where it is 
most accessible to observations. During this epoch large- 
scale (Mpc) ionized (HII) regions grow around clusters of 
sources and delineate a density contrast of order unity in 
the hydrogen abundance, far greater than the underly- 
ing fluctuations in the matter density on the same scales 
44]. Analytic models of the HII bubble size and re- 
sulting^ lem correlations were formulated by Furlanetto 
et al. [25| , who considered bubbles to self-ionize once they 
contained a sufficient mass fraction in galaxies (ionizing 
sources). This basic formalism has been applied to nu- 
merical simulations [45| for more detailed calculations of 
the ionization power spectrum. Complementary simula- 
tions by Trac and Cen [46| have been compared directly 
to the Furlanetto et al. [25j model showing reasonable 
agreement 21]. For computational ease, we choose to 
use the fully analytic calculation of 25] for the ioniza- 
tion fluctuations, including a correction for an error in 
the bubble bias as noted by McQuinn et al. 47]. This 
prescription allows us to calculate P xx and P x s- The 
key parameter in the Furlanetto et al. [25[ model is the 
ionizing efficiency of sources £; we extract this redshift- 
dependent quantity from our astrophysical model by re- 
lating the ionized fraction Xi to the total collapse fraction 
/con through Xi = C/coU- 

Fluctuations in the Lya coupling originate from spa- 
tial variation of the Lya flux, which arise primarily due 



to the strong clustering of early sources. The calcula- 
tion of these fluctuations has been discussed in detail 
by Barkana and Loeb (2(| and Pritchard and Furlanetto 
|28j |. We consider fluctuations from L ya photons arising 
from both stellar and X-ray sources [27], IH, [49[, which 
have different scale dependence. We follow the model of 
Barkana and Loeb [26[ in calculating S a by integrating 
the flux from all sources weighted appropriately by dis- 
tance and in addition account for higher Lyman series 
photons [28l . l38j . Recent investigations have shown that 
including the scattering of Lya photons in the IGM near 
to sources redistributes the Lya flux leading to greater 
power on large scales [5(1 HH, HI]. We neglect this ef- 
fect estimating that this would change our results at the 
10 — 20% level when Lya fluctuations dominate. 

Temperature fluctuations arise initially through the 
competition between adiabatic cooling and Compton 
scattering on the CMB (which drives the gas towards 
being isothermal). After thermal decoupling at z ~ 200, 
the former process dominates and temperature fluctua- 
tions grow since the rate of cooling is density dependent. 
Once star-formation begins, X-ray heating drives large 
temperature variation due to its highly inhomogeneous 
nature. We calculate the evolution of temperature fluc- 
tuations using the formalism of Pritchard and Furlanetto 
[291 ] including all three processes. Temperature fluctua- 
tions may also arise as a result of exotic heating mech- 
anisms such as the decay of dark matter [3]. In this 
paper, we are interested in producing a "standard" pic- 
ture of the 21 cm signal and will leave exotic deviations 
to future work. 

Before star-formation, 21 cm fluctuations are deter- 
mined by simple atomic physics and so, potentially, pro- 
vide a window into the physics of inflation, via precise 
measurements of the spectrum of density fluctuations 
(l2T |. Detailed numerical work [l?], [H| has shown this 
period to be rich in physics as recombination, thermal 
coupling, gas pressure and other physical processes in- 
teract to determine the final 21 cm signal. Since Comp- 
ton drag on the CMB distinguishes between baryons and 
dark matter, detailed calculation must track the evo- 
lution of both components separately [54|. After ther- 
mal decoupling, however, baryons fall into the potential 
wells created by the dark matter and by z < 30 it is a 
good approximation to take the baryons as tracking the 
dark matter, although this ignores the presence of acous- 
tic oscillations in the baryon component and pressure 
smoothing on small scales. Baryon acoustic oscillations 
(BAO) represent a major target for galaxy surveys and 
21 cm observations since they help to constrain dark en- 
ergy via the angular-diameter-distance redshift relation 
[551 ]. 21 cm observations of baryon oscillations comple- 
ment galaxy surveys well [HI, [H, [54| , since they cover a 
higher redshift range and have a different set of biases 
than galaxy surveys which may be affected by inhomo- 
geneous galaxy formation [56j . 

In our calculations, we assume that baryons trace the 
underlying density field on large scales but are smoothed 
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by the finite pressure of the gas on small scales [57j so 
that S b = 5(1 + fc 2 i?p) _1 , where R F = l/k F and k F is 
the filtering scale [58(. Thermal broadening introduces 
a second smoothing scale Rt [26j], which we incorporate 
through a Gaussian cutoff factor, exp(— k 2 R^). Until 
X-ray heating becomes important R F > Rt, while af- 
terwards Rt > R f as the filter mass takes a finite pe- 
riod of time to respond to the greatly increased gas tem- 
perature. These smoothing scales are small (typically 
Rt ~ R F ~ kpc) until photoionization heats the gas 
to ~ 10 4 K and smoothing becomes important for 
wavenumbers k > 10 Mpc -1 . We calculate linear density 
fluctuations using the fitting function of [59[ . Non-linear 
density effects become important even at high redshift 
and we incorporate these using the halo model [60j ] . 

After reionization is complete, a few percent of all neu- 
tral hydrogen (by mass) remain in over dense clumps that 
are self-shielded (CI hi = 10 -3 ). These clumps are ob- 
served as high-column density absorbers in quasar spec- 
tra in the form of Lyman-limit or damped Lya systems 
flU ]. Fluctuations in the 21 cm brightness from this post- 
rcionization epoch will trace the density fluctuations with 
a nearly constant bias fl8l|. Following the calculations 
of Wyithe and Loeb we take the post-reionization 
signal to be P 2 i(k) = b 2 21 P S5 (k) with b 2 i = 0.03 K. Al- 
though the fluctuations here are much weaker than before 
reionization, the vastly diminished foreground at the ap- 
propriate higher frequencies makes observations feasible. 
Observation of fluctuations in this regime may be used to 
detect baryon acoustic oscillations in the matter power 
spectrum and constrain the equation of state of the dark 
energy at higher redshifts than those accessible to super- 
nova or galaxy surveys [16| . 

The existence of numerous cross-terms in Eq. (|7|) con- 
necting different fluctuations requires some thought. Fol- 
lowing the approach of Barkana and Loeb [26(, we relate 
S a to the baryon density field using a scale-dependent 
linear function 8 a — W a (k)Sb. Similarly we set 5t — 
gr(k)5b [29(. In this way, we are able to relate the cross- 
terms Pqt, PbT, Pba, and equivalent 8 terms to Pss- 
Since the Lya flux, gas temperature, and baryon density 
do not interact directly with one another these cross- 
terms should be a good approximation. 

The treatment of the neutral fraction fluctuations is 
more complicated. Once HII bubbles begin to form, we 
are faced with a two phase medium composed of almost 
fully ionized HII regions and the largely neutral IGM out- 
side. We must distinguish between x it the fraction of the 
cosmic gas incorporated in fully ionized bubbles, and x e 
the ionized fraction in the largely neutral IGM (which 
could be increased beyond its primordial value due to 
diffuse X-rays). Neutral fraction fluctuations arise from 
both, although in most models that we consider x e re- 
mains small at all times and so fluctuations in x e are also 
small. The calculations of S a and St by_ Barkana and 
Loeb [26| and Pritchard and Furlanetto [29| assume that 
the filling fraction of HII regions is negligible. However, 
here we must relax that assumption in order to calculate 



terms like P xa - Since there is essentially no 21 cm sig- 
nal from the ionized HII regions, the HII regions act to 
mask the fluctuations in a and T. Overdense regions will 
host sources of Lya and X-ray photons, as well as ioniz- 
ing photons, suggesting that cross-terms of the form P xa 
and P x t may be complicated. We show later that these 
cross-terms are subdominant, so that, for a first approx- 
imation, we may neglect these complications and model 
the cross-terms as P xa = W a (k)P xb , P xT = g T (k)P x b, 
and so on. Detailed calculation will require numerical 
simulation of the ionization and radiation fields to deter- 
mine the cross-terms. 

We apply the model of Furlanetto et al. [25| to calcu- 
late P xx , P x s, and P x s x s- Note that this model first cal- 
culates the correlation functions and Fourier transforms. 
For the quartic term, we ignore the possible connected 
part so that P x g x s — P xx Pss + 2Pj ! (5 . Quartic terms in- 
volving 5 V contain angular dependence that threatens the 
possibility of the ji 2 decomposition [20| . The other quar- 
tic terms we handle in the same way as terms like P xa , 
first reducing them to a quartic of 5 and X multiplied by 
factors of W a and gx, then calculating them in the above 
way. These are generally small corrections, relevant only 
at the end of reionization. 



III. RESULTS 
A. Mean History 

We choose three astrophysical models chosen to repro- 
duce the WMAP3 central and ±1 — a for the cumula- 
tive electron-scattering optical depth r. Model A uses 
(MonjGM, fa, fx, /*) = (200,1,1,0.1) giving z rcion = 6.47 
and r = 0.063. Model B uses (Mon.iGM, fa, fx, /*) = 
(600,1,1,0.2) giving z rc ion = 9.76 and r = 0.094. Model 
C uses (AWigm, fa, fx, /*) = (3000,0.46,1,0.15) giving 
Zreion = H-76 and r = 0.115. The evolution of the mean 
temperatures with redshift for these models is shown in 
Figure [U The top axis shows the relevant observing fre- 
quency range given by v h s = v^icm/ (1 + z )- Since we are 
leaving the cosmology unchanged between models, the 
mean history only varies once star-formation becomes 
important around z « 25. The most significant differ- 
ence between these models is the redshift of reionization 
that cuts off the 21 cm signal at low redshift and is most 
responsible for the different r values. Model C has a sig- 
nificantly decreased Lya flux, which shows itself through 
a decrease in the strength of the absorption signal. 

The global evolution of the 21 cm signal is itself the 
target of several experiments, e.g. the "Cosmological 
Reionization Experiment" (CoRE) [6l[ and the "Exper- 
iment to Detect the Global EOR Signature" (EDGES) 
[62l ]. Although such an observation is conceptually sim- 
ple, it is experimentally challenging. Distinguishing be- 
tween a global signal and other sources of all sky emis- 
sion, including Galactic synchrotron, free-free radiation, 
and the CMB is difficult. Experimental detection re- 
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FIG. 1: Top panel: Evolution of the CMB temperature Tomb 
(dotted curve), the gas kinetic temperature Tk (dashed curve), 
and the spin temperature Tg (solid curve). Middle panel: 
Evolution of the gas fraction in ionized regions Xi (solid curve) 
and the ionized fraction outside these regions (due to diffuse 
X-rays) a;,; (dotted curve). Bottom panel: Evolution of mean 
21 cm brightness temperature TJ,. In each panel we plot curves 
for model A (thin curves), model B (medium curves), and 
model C (thick curves). 



lies upon reionization proceeding rapidly leading to a 
distinctive step-like feature in the frequency direction, 
which would not be expected to be produced by the 
spectrally-smooth foregrounds. With the assumption of 
sharp reionization, EDGES [62j places an initial con- 
straint that Tb < 450 mK at z = 8. While this is far from 
the expected signal amplitude, such constraints will im- 
prove with time. Efforts are also underway to extend the 
frequency coverage to v w 50 MHz to access the transi- 
tion from an absorption to emission signal. As Figure Q] 
indicates, this transition is likely to be significantly larger 
in amplitude (~ 100 mK) than that at the end of reion- 
ization (~ 20 mK). 



B. Fluctuation History 

The three dimensional nature of the 21 cm signal 
makes it difficult to convey the evolution of the fluctua- 
tions with a single 2-dimensional plot. We therefore plot 
the evolution of four individual comoving wavenumbers 
k = 0.01, 0.1, 1, and 10 Mpc -1 , spanning the range 
that might be observed. On large scales we expect con- 
tamination from foregrounds to limit the detection of the 
power spectrum. On small scales thermal broadening of 
the 21 cm line will smooth the signal. It is also to be ex- 



pected that many of our approximations will break down 
as small scale information about the sources becomes im- 
portant (see for example 63] for the importance of higher 
order correlations on small scales during reionization). 
For the mean histories shown in Figure [1] we calculate 
the evolution of the 21 cm angle-averaged power spec- 
trum, which is plotted in Figures fJl EH and |H for models 
A, B, and C, respectively. 

The evolution of At 6 clearly shows three regimes: the 
post-reionization regime at low redshifts (z < z ro ion) 
where the 21cm fluctuations from residual hydrogen fol- 
low the matter power spectrum, an intermediate red- 
shift regime (x re ion < z < ztrans) where Lya coupling 
produces a large signal and complicated astrophysics 
leads to significant scale dependence, and a high red- 
shift collisionally-coupled regime where 21 cm fluctua- 
tions track the density field (z > z trans « 23). For peda- 
gogical purposes, let us describe the evolution on a single 
comoving scale (say, k = 0.1 Mpc -1 ) and draw attention 
to the main features. Thermal decoupling at z ~ 200 is 
a gradual process and, initially, Ax b grows due to a com- 
bination of the growth of density fluctuations and the 
steady gas cooling below 1cmb- As the gas rarifies and 
cools, collisional coupling becomes less effective and, at 
z ~ 60, Ax b begins to decrease in amplitude. Note that 
the continuing growth of structure offsets the turnover 
on A.T b from the minimum of seen in Figure [1] to oc- 
cur at z ~ 90. As collisional coupling diminishes, the 
signal drops towards zero. This occurs while Tk < 30, 
a regime where k\-q{Tk) drops exponentially with Tk 
[12] and results in a rapid drop of the signal at z < 40. 
Before the signal drops all the way to zero, significant 
star- format ion occurs and the resultant Lya production 
leads to the beginning of Lya coupling by z rs 25. The 
exponential increase in the global star formation rate at 
these redshifts is responsible for the rapid increase in T\, 
and Ax b - With this rise in signal, we enter into a regime 
dominated by astrophysics and begin to see complicated 
scale dependence. 

Initially, Lya fluctuations boost the signal somewhat 
above the level of density fluctuations alone. However, X- 
ray heating follows not far behind and contributes to Ax b 
with the opposite sign (hotter regions produce a weaker 
absorption signal, see Pritchard and Furlanetto In 
this competition, X-ray driven temperature fluctuations 
dominate causing Ax b to pass through a zero point (seen 
as a sharp dip at z ~ 18 in all three plots). Tempera- 
ture fluctuations dominate as Tk approaches Tomb and 
Tb vanishes. In proceeding to the emission regime, we 
note based on Figure [1] that the brightness fluctuations 
emitted Tb are generically smaller than they were during 
the absorption regime, leading to a decreasing trend in 
A.T b . As reionization gets underway, ionization initially 
causes Ax b to drop leading to a pronounced dip in its 
evolution. This occurs as a result of the clustering of 
ionizing sources in over dense regions causing the ionized 
HII regions to "mask out" those dense regions that have 
been producing the strongest 21 cm signal. As reioniza- 



7 



tion proceeds, the contrast between ionized and neutral 
regions comes to dominate and Ay b , rises until xh ~ 0.5 
after which the contrast begins to drop. 

Towards the end of reionization the signal drops 
sharply as very little gas is left neutral. The post- 
reionization signal grows slowly as the density field grows. 
Since by this time the gas is photo-heated to Tjc ~ 
30, 000 K the thermal width of the 21 cm line is sufficient 
to smooth out the signal on wavenumbers k > 10 Mpc -1 . 
This cutoff potentially acts as a thermometer of the gas 
after reionization giving information about the tempera- 
ture of gas contained in dense clumps. 

As a result of the interplay between the radiation 
fields, as A.T b evolves it shows three peaks within the 
astrophysics-dominated regime. An important feature of 
this complicated evolution is that the maximum ampli- 
tude of Ax,, occurs at different k values for different red- 
shifts. Accurate observation and modeling of this com- 
plicated evolution may provide important information 
about the early radiation fields. 
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FIG. 2: Redshift evolution of the angle- averaged 21 cm power 
spectrum At 6 for Model A at k — 0.01 (solid curve), 0.1 (dot- 
ted curve), 1.0 (short dashed curve), and 10.0 (long dashed 
curve) Mpc -1 . Reionization at z = 6.5. 
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FIG. 3: Redshift evolution of the angle- averaged 21 cm power 
spectrum Ar s for Model B. Reionization at 2 = 9.8. Same 
line conventions as Figure [5] 
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It is helpful to get a sense of how the amplitude of the 
signal compares with galactic foregrounds. We take the 
sky noise to be T sky « 180 K(^/180MHz) -2 - 6 (appropri- 
ate for galactic synchrotron emission [1(|), noting that 
the normalization depends upon the region of sky being 
surveyed. In Figure [21 Figure [31 and Figure [H we plot 
rT s ky(f) where r ranges from 10 -4 — 10 -9 . We see that 
reducing foregrounds by a factor of ~ 10 -5 is required to 
observe fluctuations during reionization and cosmic twi- 
light. The difficulty increases if reionization occurs early, 



FIG. 4: Redshift evolution of the angle- averaged 21 cm power 
spectrum At 6 for Model C. Reionization at z = 11.8. Same 
line conventions as Figure [2] 



which has the effect of compressing the signal at high red- 
shifts (model C). The signal from astrophysics in these 
three models begins at v w 60 MHz and continues to 
v w 150 MHz although this upper limit is very sensitive 
to the details of reionization. 

Removing foregrounds to the rather low level of ~ 10 -7 
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is required to have a hope of observing the Dark Ages in 
the range z = 30-50 (25 MHz < v < 60 MHz). Some 
hope is restored by noticing that the rapid evolution of 
Ax b between z = 30 — 50 mimics the strong frequency 
dependence of the foregrounds. If one can remove fore- 
grounds to the level needed to observe z = 30, then the 
same level of foreground subtraction would allow z = 50 
to be observed. Beyond z = 50, the amplitude of Ar b 
decreases rapidly and becomes progressively more diffi- 
cult to observe. Modifying the thermal history between 
z = 30 and recombination by introducing exotic physics 
can greatly change this high redshift behavior. Energy 
injection by decaying dark matter [14j can increase the 
residual electron fraction in the IGM while maintaining 
collisional coupling to lower redshifts boosting the high 
redshift 21 cm signal and improving the chance of ob- 
serving the "Dark Ages". 



in the middle panel of Figure [5] to illustrate where dif- 
ferent effects become important. Since these terms dif- 
fer only in one element of the cross-correlation, A M 2 is 
more sensitive to the form of the different fluctuations 
than A„o. At z > 22, the dominant contribution to A^2 
comes from the bV term although the TV term is also 
significant. The TV term changes sign at z w 38 owing to 
the behavior of /3y. Below z w 22 astrophysics becomes 
important and the aV term dominates for a short while 
before the TV term overtakes it, causing A M 2 to change 
sign. Around z « 15 all components drop towards zero 
briefly as Tj, = 0. This imprints an extra dip in the A^2 
power spectrum that is not seen in A M o, where the TT 
term dominates at this point. Utilizing the different red- 
shift evolution of these power spectra is likely to greatly 
increase the chance of separating different astrophysical 
sources of fluctuations. 



C. Decomposition of fluctuations 

The total 21 cm signal is composed of fluctuations in 
five quantities. Thus, when we form the power spectrum 
from Equation ([7]) we must calculate 15 distinct correla- 
tions and cross-correlations. In addition, as fluctuations 
in the neutral fraction can be of order unity, there are a 
number of fourth order correlations which must be calcu- 
lated. Separating out these different components will be 
difficult. A first step will be to exploit the angular depen- 
dence of the peculiar velocity term to separate out the 
contributions from P„o , P^i , and P^i . Since each of these 
three terms depends on different combinations of terms, 
measuring the three separately gives extra discriminating 
power [23| . This decomposition is complicated by the 
presence of the P X xS b s v and P X p5„&., terms, which show 
very different /i dependence [20]. Since these terms will 
be large only once reionization is well underway, we shall 
ignore the angular dependence of these terms. 

We plot the different angular power spectra for model 
A in Figure [51 The fj, 4 term is the simplest, arising solely 
from the VV term, so that it tracks Pss, making it ideal 
for trying to measure cosmology. Its evolution is modu- 
lated by the evolution of the mean Tj,, being divided into 
two absorption bumps and an emission bump. 

The /Lt contribution contains contributions from 10 
different terms, and so we have not attempted to show 
its full decomposition. Instead, we plot the contribution 
only from the auto-correlations, which gives some idea of 
where different contributions are important. Note that 
the contribution of cross-terms is significant and that 
many cross-terms have negative signs over some range 
of redshifts. Hence, A„o can be smaller than the sum of 
the auto-correlations would indicate. Since the /jP term 
dominates Ax b , its behavior is very similar in form to 
that plotted in Figure [U Temperature, Lya , and ioniza- 
tion show clear overlapping regions of contribution. 

Finally, the /i 2 contains contributions from the bV, 
XV, aV, and TV cross-terms. We have plotted these 
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FIG. 5: Angular decomposition of 21 cm power spectrum for 
model A at wavemmiber k = 0.1 Mpc - . Top panel: A (J o(fe) 
(solid curve). We also show the separate contribution from 
the bb (dotted curve), aa (long dashed curve), TT (short 
dashed curve), and XX (dot-dashed curve) terms. Note that 
there are additional cross-terms, many of which contribute 
to the power spectrum with negative sign, not shown here. 
Middle panel: A M 2(fc) (solid curve). We also plot the contri- 
bution from the bV (dotted curve), aV (long dashed curve), 
TV (short dashed curve), and XV (dot-dashed curve) terms, 
indicating the sign of the contribution to the power spectrum 
as positive (thick curves) or negative (thin curves). Bottom 
panel: A 4 (A;) (solid curve). 
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D. Separation of fundamental physics from 
astrophysics at high redshifts 

We next turn our attention to characterizing the red- 
shift at which the 21 cm signal becomes affected by the 
astrophysics of star formation. Although, in principle, 
this could occur when any of /3tSt, ftx&x, or (3 a S a be- 
come large in comparison with faSb, as we have seen in 
the preceding section, the onset of Lya coupling is likely 
to provide the first astrophysical modification of the sig- 
nal. In order to quantify this transition, we will con- 
sider the redshift ztrans.a at which (3 a 5 a = f3t5b for the 
first time. Although cosmological parameters may be ex- 
tracted from components other than the 66-term, this is 
the dominant component at high redshift. 




0.01 0.1 1 10 



f«. fx 

FIG. 6: Evolution of ztrans.a with f a (thin curves) and fx 
(thick curves) for k = 0.01 (dotted curves), 0.1 (short dashed 
curve), and 1 Mpc -1 (long dashed curve). Also plotted is 
the redshift at which x a — x c (solid curves). 



Figure [6] shows the dependence of ztrans,a on the pa- 
rameters fx and f a . In order to separate out the differ- 
ent effects, we consider two artificial situations. First, we 
vary /„ using the fiducial value fx = 1, but neglecting 
Lya photons produced by X-rays. Next, we vary fx and 
neglect the stellar Lya component, i.e. f a = 0. Thus, 
the plotted values of Ztrans.a should be regarded as lower 
limits, since including the Lya flux of the missing com- 
ponent will cause the transition to occur slightly earlier. 
All other parameters are set to those of model A. We 
plot ztrans.a for three different wavenumbers, k = 0.01, 
0.1, and 1 Mpc -1 , noting that there is little difference 
between the three curves. For comparison, we also plot 
the redshift at which Lya and collisional coupling are 
comparable, defined by x a = x c . 



The basic conclusion is that for a Gaussian field of 
initial density perturbations, ztrans.a depends only log- 
arithmically on f x and f a . This dependence arises be- 
cause the fraction of mass in galaxies grows exponen- 
tially with cosmic time on the high-density tail of the 
Gaussian. Exploring four orders of magnitude in f a 
causes ztrans.a to vary by a redshift interval of Az — 8 
(A^ ks 20 MHz). The same variation in fx causes ^trans.a 
to vary by Az = 14 (Av w 40 MHz). This extra depen- 
dence arises because modifying f x also changes T K . Note 
that we can relate the uncertainty in the star formation 
efficiency /» to these results, provided we remember that 
changing /* will modify the production of both UV and 
X-ray photons. 

We can understand this parameter dependence, by 
considering the redshift at which Lya coupling first be- 
comes important. Obtaining x a = 1 requires a flux 
J a = 0.0767 (1 -I- z/20) 2 Lya photons per baryon. We 
can convert this into the fraction of baryons in stars 
needed for Lya coupling 

To a good approximation, we can find the redshift of 
transition by setting rj a = f co \\f* to find the point when 
enough baryons are in stars for the transition to occur. 
This will slightly overestimate the redshift of transition, 
as we are not accounting for photons that redshift out of 
the Lya resonance. Since / co u initially grows exponen- 
tially with redshift we obtain a logarithmic dependence 
on f*N a . A similar result will apply for Lya photons 
from X-ray sources. 

In more extreme astrophysical models, it is conceivable 
that intense X-ray heating might affect the 21 cm signal 
before Lya coupling became important. We therefore 
also consider the redshift Ztrans.T at which /3tSt — Pb$b 
for the first time. This is plotted in Figure[7]for the same 
models used above. For comparison, we also plot the red- 
shift at which the gas is heated to Tk — Tcmbj which 
occurs significantly later and marks the transition from 
absorption to emission. Unsurprisingly, z tranSj T shows 
very little dependence on f a . The dependence on fx is 
more pronounced. In no part of the parameter space ex- 
plored does ztrans,T exceed Ztrans.a, suggesting that Lya 
coupling will always be the first source of astrophysical 
fluctuations. This need not be the case if gas heating 
does not also produce Lya photons, for ex amp le if shock 
heating dominates in the early Universe [6J|. We also 
note that values of f*fx ^ 0.01 are required to ensure 
that Tjc > Tcmb m the redshift range z < 12 that will 
be probed by the first low-frequency observatories. 

In the same fashion as the onset of Lya coupling, we 
can seek to understand the requirement for heating the 
gas above the CMB. Assuming that Tk <C Tomb at the 
onset of heating, an energy of AE w 3fcsTcMB/2 per 
baryon is required. Using parameters appropriate for X- 
ray emission for starburst galaxies and assuming that x e 
is approximately the primordial value, yields the fraction 
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FIG. 7: Evolution of ztrans,T with the same line conventions 
as Figure [6] Also plotted is the redshift at which Tk = Icmb 
(solid curves). 



IV. OBSERVATIONAL PROSPECTS 

In section fill B[ we showed that detecting the 21 cm 
signal would require foreground removal to the level of 
10 -5 for reionization and cosmic twilight and 10~ 7 to 
reach the Dark Ages. Next we explore the experimen- 
tal capabilities necessary to achieve the required sensi- 
tivity for an accurate measurement of the 21cm power 
spectrum. Before presenting our calculations, let us di- 
gress and consider the angular scales probed by these 
experiments. Until this point our calculations have been 
presented in terms of comoving wavenumber k, appro- 
priate for a 3D power spectrum. For observers the 
more natural quantity is the angular scale on the sky, 
0. Figure [5] shows a conversion between the two using 
= (2n / k) / cIa(z) , with cIa(z) the angular diameter dis- 
tance to redshift z. We also relate angles to multipole 
index t in the spherical harmonics decomposition of the 
sky, based on the approximate relation I w 7r/0. Al- 
though this conversion is somewhat crude, it allows a 
quick comparison with quantities that might be more in- 
tuitive to the CMB researchers. Astrophysics primarily 
affects the 21 cm power spectrum on comoving scales 
k = 0.1 — 1 Mpc -1 corresponding to angular scales of 
several arcminutes. 



of baryons in stars required to heat the gas 
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(ii) 

Comparison with Eq. (JTUJ) shows that heating the IGM 
to Tcmb with X-rays requires considerably more star for- 
mation than producing Lya coupling. While generalizing 
this argument rigorously from the mean history to the 
case of fluctuations is difficult, the basic point remains. 
Finally, in the absence of significant clumping reioniza- 
tion requires approximately one UV photon per baryon 
in stars so that we may write 
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= 2.5 x 10" 



ion. IGM 



/4000\ /0.1 \ 



(12) 



Thus reionization should occur after Lya coupling is sig- 
nificant and the gas has been significantly heated. 

We conclude that for astrophysical parameters not too 
far from our fiducial ones fx ~ fa ~ 1 and /* m 0.1, 
Lya fluctuations will first become important at z ss 25 
(v w 55 MHz) and temperature fluctuations at z ~ 19.5, 
although this latter number is considerably more uncer- 
tain. This defines the minimum requirement for a low- 
frequency observatory hoping to probe the cosmology 
dominated regime before the first sources affect the 21 
cm signal. For the mean history, Lya coupling becomes 
relevant at z rj 23 and Tk > Tcmb at z < 14. These 
results are somewhat sensitive to the cosmology used, 
especially the value of er 8 , which affects the collapse frac- 
tion. The location of these landmarks will be important 
for attempts to probe the mean 21 cm signal. 
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FIG. 8: Redshift evolution of angular scales. B = 
(2n / k) / cLa{z) , which we map to angular moments using I « 
7r/0. This is meant to be an approximate guide, rather than 
an exact conversion. 



We now calculate the signal-to-noise ratio (S/N = 
\J -P21 /SP21 ) for three fiducial experiments: (i) a 
pathfinder class experiment like the Murchison Widefield 
Array (MWA); (ii) a fully fledged EoR instrument like 
the Square Kilometer Array (SKA); and (Hi) a futuristic 
lunar array (LA) (65| . The parameters used to describe 
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these instruments are summarized in Table Q] In each 
case, we specify the number of antennae N a , and the 
total collecting area j4tot — N a A e for an experiment op- 
timized to observe z = 8. We also specify the bandwidth 
B and total observing time tint ■ Note that while we have 
chosen parameters that correspond roughly to proposed 
experiments, we take liberties with the antennae config- 
urations to maximize the sensitivity of each experiment 
at each redshift. We are interested in the most favorable 
capabilities of each experiment given the parameters in 
Table 1 , and so we optimize the associated array configu- 
ration for each redshift window. For ease of reference, we 
will label these experiments MWA, SKA, and LA. Cur- 
rent designs for the actual MWA and SKA limit their 
frequency range (> 80 MHz) and antennae distribution, 
and so the above labels should not be associated with the 
actual instruments being designed or constructed. The 
labels are simply meant to denote different scales of ex- 
perimental effort. 

TABLE I: Low-frequency radio telescopes and their parame- 
ters. We specify the number of antennae N a , total collecting 
area A to t, bandwidth B, and total integration time tint for 
each instrument. 



Array 


N a 


At oi (W°n?) 


B (MHz) 


tint (hr) 


MWA 


500 


7.0 


6 


1000 


SKA 


5000 


600 


6 


1000 


LA 


7800 


3600 


8 


12000 



The variance of a 21 cm power spectrum estimate for 
a single k-mode with line of sight component k\\ — \xk is 
given by [6?| : 

—— T 2 P 21 (k^)+T 2 ys — nTY ( ~~a~ ) ■ 

Aficid [ y Bt int n{k±) \A e J 

(13) 

We restrict our attention to modes in the upper-half 
plane of the wavevector k, and include both sample 
variance and thermal detector noise assuming Gaussian 
statistics. The first term on the right-hand-side of the 
above expression provides the contribution from sample 
variance, while the second describes the thermal noise of 
the radio telescope. The thermal noise depends upon the 
system temperature T sys , the survey bandwidth B, the 
total observing time t- m ti the conformal distance D(z) to 
the center of the survey at redshift z, the depth of the 
survey AD, the observed wavelength A, and the effective 
collecting area of each antennae tile A e . The effect of 
the configuration of the antennae is encoded in the num- 
ber density of baselines n± (k) that observe a mode with 
transverse wavenumber k± [20j |. Observing a number of 
fields Afield further reduces the variance. 

Estimates of the error on a power spectrum measure- 
ment are calculated using the Fisher matrix formalism, 



so that the 1 — a errors on the model parameter Xi are 
(Fr. 1 ) 1 /2 ; wri ere 

= \ - ek 3 V smvcy 1 dP Tb dP Tb , . 

y ^ 4^ 2 <J 2 P (k,n) dX l d\j ' 1 ' 

In this equation, Purvey = D 2 AD(X 2 /A e ) denotes the 
effective survey volume of our radio telescopes and we 
assume wavenumber bins of width Ak = efc. We will 
be interested in the cases where Xi — {Pr b } and Xi — 
{P^,P^,P^}. 

Instrumental sensitivities have been shown to depend 
upon the distribution of antennae. One typical model is 
to assume a close packed core with filling fraction close 
to unity surrounded by antennae distributed with a r~ 2 
dependence out to a sharp cutoff at the edge of the array. 
Lidz et al. [6q | showed that significant gains in sensitiv- 
ity can be achieved by using a "super core" configuration, 
condensing the array so that the filling fraction is unity 
throughout the array. Since most of the longest baselines 
in the r~ 2 configuration only poorly sample the largest 
k modes anyway, this enhances the signal-to-noise while 
losing measurements on only a few large k modes. Since 
we wish to maximize the sensitivity of our instrument, 
we will use a similar "super core" configuration. The ef- 
fective collecting area of the antennae is a strong function 
of observed frequency (scaling as v~ 2 or A 2 ), so long as 
antennae are more widely spaced than A/2. We optimize 
our arrays at each redshift by setting the minimum base- 
line to A/2 and then close packing the antennae. Our 
calculations should then give a reasonable idea of the 
best performance that can be hoped for given our chosen 
specifications. In practice, an array will be optimized for 
a single redshift and will suffer performance degradation 
when probing higher redshifts, where geometrical shad- 
owing becomes important. Additionally, some fraction 
of antennae will be reserved to provide the long baselines 
needed for foreground removal. 

Figure [5] shows three redshift slices of the 21 cm power 
spectrum for model A. These illustrate the shape of the 
sensitivity of the various instruments. The sensitivity 
is degraded by foregrounds on large scales and thermal 
noise on small scales. In each case, those scales closest to 
the foreground cutoff (which is determined by the line- 
of-sight resolution set by the experimental bandwidth) 
on large scales are the most likely to retain some fore- 
ground contamination. Unfortunately, these tend to be 
the modes that are measured with the highest sensitivity. 
The shape of the power spectrum on intermediate scales 
shows significant evolution as different sources of fluctu- 
ations become important. At z = 30.2, Ax b is dominated 
by density fluctuations. At z = 15.7, temperature fluc- 
tuations lead to a trough on a scale k » 1 Mpc -1 and 
a peak at k « 0.1 Mpc~L The location of this trough 
evolves at higher redshift to lower k, where it might be 
more easily detectable. By z = 8, ionization fluctuations 
flatten the power spectrum significantly. The evolution 
of the shape of the power spectrum during reionization 
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has been examined in more detail by Lidz et al. [66| , who 
explored the ability of the MWA to measure the slope 
and amplitude of the power spectrum. While the SKA 
is a significant improvement over the MWA, it takes the 
futuristic LA to directly observe the trough in Ax b at 
z = 15.7 or to detect fluctuations at z — 30.2. 
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FIG. 10: Signal to noise ratio as a function of redshift for 
our optimized experiments. Curves indicate k — 0.1 (dotted 
curve), 1.0 (short dashed curve), and 10 Mpc -1 (long dashed 
curve). Top panel: MWA. Middle panel: SKA. Bottom 
panel: LA. 



FIG. 9: Errors on the angle-averaged power spectrum At 6 
(solid curves) for our optimized experiments MWA (dotted 
curve), SKA (short dashed curve), and LA (long dashed 
curve), as functions of wavenumber k at redshifts z = 7.9 
(top panel), 15.7 (middle panel), and 30.2 (bottom panel). 



Figure [TU] shows achievable S/N ratio for measure- 
ments of the spherically averaged signal Ax b for the 
three instruments as a function of redshift for three 
different wavcnumbers k — 0.1, 1, and 10 Mpc -1 . 
These wavenumbers span the scales to which these in- 
struments are sensitive. We see that a well optimized 
MWA achieves high signal-to-noise ratio S/N for only 
the k = 0.1 Mpc -1 mode and moderate S/N for k 
1.0 Mpc -1 . SKA achieves the sample variance limit at 
k = 0.1 Mpc -1 over a wide range of scales and shows sen- 
sitivity to k = 1.0 Mpc -1 through most of the EoR. Only 
the LA is capable of detecting the signal from the Dark 
Ages at z > 30. These curves are optimistic since the 
instruments are assumed to be optimized to the specific 
frequency of observation. At very low redshifts (z < 3), 
the bandwidth that we have chosen is inadequate to the 
task of removing foregrounds on the scale k = 0.1 Mpc -1 , 
hence the sharp cutoff in sensitivity seen in Figure [TU1 

We now turn to the possibility of performing the angu- 
lar separation by powers of /j, discussed earlier. There are 
two main reasons for performing this separation. First, 
as discussed in section ^IIII the 21 cm fluctuations aris- 
ing from different astrophysical models tend to overlap, 
giving Ax b a complicated shape. In order to constrain 



astrophysics, it may be advantageous to turn to the A^2 
term. This term shows more structure than the spheri- 
cally averaged signal (see Figure [TTj) . while being simpler 
to analyze as it is composed of fewer terms. Second, 
the complexity of the astrophysical fluctuations obscures 
the underlying cosmological information contained in the 
density contribution to 21 cm fluctuations. This problem 
can be circumvented by measuring A M 4, which is deter- 
mined by the density field alone [231] . Although the use 
of this term is most important during the astrophysics 
dominated regime, it might also be important during the 
Dark Ages. If exotic energy injection mechanisms, such 
as decaying dark matter [14J , occur during the Dark Ages 
they would produce fluctuations in temperature and ion- 
ization fraction that would interfere with the extraction 
of cosmology from the density fluctuations. The angular 
separation would provide extra information about these 
processes and still allow cosmological constraints to be 
extracted. 

Figure [TT] shows the sensitivity of our fiducial experi- 
ments to the angular components as a function of redshift 
for wavenumber k — 0.1 Mpc -1 . All three experiments 
show some sensitivity to A M o , which is to be expected 
from our analysis of the spherically averaged signal. Only 
SKA and LA show any sensitivity to A„2 and A^4 . 

Since A M o dominates the signal, the errors in A M 2 and 
A^4 track its shape. This results in only a few windows 
where observation of these components is possible. Dur- 
ing the end of reionization the largely isotropic ionization 
signal dominates [201 ] making observations of the /x 2 and 



13 



6 



6 



10 


- k = 0.1 Mpc" 1 f)x 


1 'I''; 


1 


rxy / V \ 

. / ,?/ .-•'*' -—7 




0.1 




-= 


10 

1 


, m 


I 1 1 1 : 


0.1 




- 


10 


i | i i i i | i i 


I 1 : 


1 


— / *(* 

/ W v V ''''' 

I V V-'''S^ — 




0.1 


r, TV^/, . i . . 


i , , , ,"l 



5 10 50 100 

1+z 



FIG. 11: Experimental sensitivity to the separation of powers. 
Plotted is the redshift evolution of the power spectrum (thick 
solid curve) at k = 0.1 Mpc - and the corresponding sensi- 
tivity for MWA (thin solid curve), SKA (thin dotted curve), 
and LA (thin dashed curve). We also plot sensitivity for LA 
with its observing time split between 16 separate fields (thin 
dot-dashed curve). Top panel: A^o. Middle panel: A M 2. 
Bottom panel: A M 4. Calculations are for model A. 



ment will gain sensitivity on scales to which it is sample 
variance limited while losing sensitivity where thermal 
noise is important. Only if the experiment has sufficient 
instantaneous sensitivity, will it be worth observing mul- 
tiple fields. Note though that for some interferometer 
designs, e.g. LOFAR, multiple fields can be observed 
simultaneously,improving the S/N while adding only to 
the computational cost of the required correlations. For 
an experiment like LA, which has the sensitivity to reach 
sample- variance limited measurements, experimental de- 
sign will be key in distributing antennae and observing 
time to minimize errors. 

Since performing the angular separation may be vital 
to obtaining both astrophysical and cosmological infor- 
mation, we briefly consider a modified version of the LA 
in which we redistribute the same fixed total observing 
time from deep integration of a single field, necessary to 
probe the Dark Ages, into shallower integrations of 16 
separate fields. This is plotted in Figure [11] (dot-dashed 
curve). We see that considerable gains in sensitivity are 
achievable, dramatically increasing the range of redshifts 
over which the angular separation is possible at high S/N 
on this scale. However, this version of the LA is unable 
to probe the Dark Ages and is sample variance limited 
over a smaller range of wavenumbers. Since observation 
strategies for direct observation of the Dark Ages versus 
separating A^4 are orthogonal to one another, detailed 
consideration of which approach is better for cosmology 
will be important for the design of future instruments. 



/i 4 components difficult. This is also true earlier when 
temperature and Lya fluctuations dominate the signal. 
Observational prospects for the angular separation are 
best when these fluctuations are small, at the beginning 
of the emission period but before reionization is well ad- 
vanced (8 < z < 12 for model A). There is also a pos- 
sibility for measuring A^2 during the absorption period 
(15<z<21 for model A) although this relies somewhat 
upon the scale-dependent cancellation between temper- 
ature and Lya fluctuations. We note that a high S/N 
detection is required for the angular separation to be fea- 
sible. For this reason, none of the instruments are able 
to perform the angular separation during the Dark Ages. 
The angular decomposition also becomes more difficult 
on small scales. 

Since the SKA and LA experiments are essentially 
sample variance limited at k = 0.1 Mpc -1 , it would seem 
that performing the full angular separation will be ex- 
ceedingly difficult except in a few redshift windows. This 
is a direct consequence of the isotropic part of the power 
spectrum dominating the uncertainty in measurement of 
the /x 2 and /z 4 terms. To improve on the results illus- 
trated here it would be necessary to beat down cosmic- 
variance by observing multiple fields on the sky. Since 

— 1/2 

the variance scales as A r ficl j while the contribution from 
thermal noise scales as \ , there is a trade-off involved in 
splitting observing time into different fields. The experi- 



V. DISCUSSION AND CONCLUSIONS 

Future low-frequency radio observations of the red- 
shifted 21 cm line hold the potential to greatly increase 
our knowledge of the high redshift Universe. Theoreti- 
cal understanding of fluctuations in the 21 cm brightness 
temperature has reached the point where it is possible to 
begin making predictions for the evolution of the signal 
from recombination to the present day. In this paper, we 
have combined the present understanding of fluctuations 
in density, neutral fraction, Lya flux, and X-ray heating, 
to calculate the 3-dimensional 21 cm power spectrum as 
a function of redshift. This draws attention to the three 
different epochs mentioned in the Introduction. Early 
reionization compresses the reionization and twilight sig- 
nal into a narrow plateau, likely complicating attempts 
to separate different fluctuations and learn about astro- 
physics. 

Our comprehensive analysis provides a useful founda- 
tion for optimizing the design of future arrays whose goal 
is to separate the particle physics from the astrophysics, 
either by probing the peculiar velocity distortion of the 
21 cm power spectrum during the epoch of reionization, 
or by extending the 21 cm horizon to z > 20 or z < 6. 

We have shown that the signal above z = 25 is likely to 
be relatively uncontaminated by astrophysics. This sets 
the minimum redshift that a low-frequency radio observa- 
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tory must probe in order to seek cosmology most simply. 
This same calculation indicates that experiments with 
sensitivity down to 60 MHz should be able to observe 
the full span of the EoR signal. 

Observations should be sensitive to both the emission 
and absorption regimes of the 21 cm signal during the 
EoR. This is in part due to the greater amplitude of 
lb during the absorption regime and in part to the en- 
hancement of the signal from temperature and Lya fluc- 
tuations. Additionally, experiments sensitive to z < 5 
should have a chance of detecting the post-reionization 
signal. 

In order for experiments to detect 21 cm fluctuations 
from the EoR they will need to first remove foregrounds 
to the level of 10~ 5 . To probe the Dark Ages, heroic ef- 
forts to obtain the level of ICE 7 will be required. If this 
can be achieved then the SKA should be able to mea- 
sure the EoR signal over several decades of wavenumber. 
Significant improvement in cither collecting area or, for 
example, going to the moon for longer integration times 
will be necessary to access the Dark Ages. 

Interestingly, the 21 cm emission from residual hydro- 
gen after reionization (z < 6) offers excellent prospects 
for probing fundamental physics [l8[, because of the 
steep decline in the Galactic foreground brightness with 
decreasing wavelength (ex A 2 6 ). The pockets of self- 
shielded hydrogen are expected to trace the distribution 
of matter, and allow a precise determination of the mat- 
ter power spectrum. This in turn will provide an un- 
precedented probe of non-Gaussianity and running of the 
spectral index of the power-spectrum of primordial den- 



sity fluctuations from inflation. Detection of small-scale 
fluctuations can also be used to infer the existence of mas- 
sive neutrinos [l3Tj and other sub-dominant components 
in addition to the commonly inferred cold dark matter 
particles. The scale of the baryonic acoustic oscillations 
in the power spectrum can be used as a standard ruler to 
constrain the equation of state of the dark energy through 
99% of the cosmic history [H, GJ, [H . 

The use of 21cm cosmology to study fundamental 
physics is not restricted to either high redshifts (z > 25; 
v < 50MHz) or low redshifts (z < 6; v > 200MHz). The 
implications for inflation and dark matter can be sep- 
arated from the astrophysical effects of star formation 
through an angular decomposition of the 21cm power 
spectrum, even during reionization (6 < z < 25). The 
angular term A^4 isolates the contribution of gravita- 
tionally induced peculiar velocities. We have found that 
with an appropriately designed instrument, separating 
the A^4 term should be possible over a wide range of 
redshifts. 
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